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bone in the framework of transversely 
isotropic poroelasticity 
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Abstract 

In this paper, a cylindrical thick-walled hollow bone filled with marrow is modeled as a transversely isotropic poroelastic 
thick-walled hollow cylinder filled with a fluid, and analytical solutions for radial vibrations in the same have been 
obtained. Considering the boundaries to be stress free, frequency equations are obtained in the cases of bone without 
marrow and bone with marrow for permeable and impermeable boundaries. Limiting the case when the ratio of 
thickness to inner radius is very small is investigated numerically. The accessible information of bony material has 
been used for the numerical evaluation. Attenuation versus ratio of outer and inner radii of bone at various anisotropic 
ratios are computed. Phase velocity and group velocity are computed as a function of wavenumber. These values are 


depicted graphically and then discussed. 
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|. Introduction 


The study of wave propagation in poroelastic solids, in 
general, is of practical importance in the fields of engin- 
eering, medicine, and geophysics, particularly, for 
investigations in destructive research areas. 
Poroelasticity has been widely used to model biological 
tissues, such as bone, cartilage, arterial walls, and brain 
as almost all tissues have some kind of fluid in their 
pores. Musculoskeltal complications, namely osteopor- 
osis (OP) not only affects bones mineral density (BMD) 
but also quality. The major consequence of OP is a 
reduction in bone density leading to bone fracture. 
The strength of bone is an important parameter for 
bone quality. Wave velocity and attenuation through 
a bone depends on the poroelastic parameters of bone, 
which are directly related to the strength of the bone. 
Bone tends to exhibit linearly elastic behavior over 
small strains, if not, it exhibits anisotropy (Humphrey 
and DeLange, 2007). Biot’s equations (Biot 1955, 1956) 
are employed to investigate the bone and the equations 
are derived from the consolidation theory to investigate 
wave propagation in the bones under the assumption 
that bone exhibits transversely isotropic behavior 


(Nowinski and Davis, 1971). In Nowinski and Davis 
(1972), bones are modeled as anisotropic poroelastic 
bodies. Cowin (1999) described fluid flow in bone tis- 
sues employing the poroelasticity theory. Analytic solu- 
tions for electromechanical wave propagation in a 
cylindrical poroelastic long bone with cavity are stu- 
died (Ahmed and Abd-Alla, 2002). Fluid flowing 
through anisotropic, poroelastic bone models in the 
opposite direction is investigated (Steck et al., 2003). 
Dispersion of axially symmetric waves in cylindrical 
bone filled with marrow is investigated under the 
assumption that bone exhibits isotropic poroelastic 
behavior (Malla Reddy et al., 2011). In Malla Reddy 
et al. (2011), a comparative dispersive study is made 
between the bone without marrow and the bone with 
marrow. Gilbert et al. (2012) made a quantitative 
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ultrasound model of the bone with blood as the inter- 
stitial pore fluid where they considered the Biot model 
treats the medium as an elastic frame with interstitial 
pore fluid. From the classical poroelasticity perspective, 
the following are a few papers on circular cylinders. 
Approximate formulas for wave lengths associated 
with the free vibrations of a thick-walled hollow elastic 
cylinder are examined by McFadden (1954). The ana- 
lytical solutions for pore pressure and stress fields for 
an inclined bore hole and the cylinder are investigated 
(Cui et al., 1997; Abousleiman and Cui, 1998). Axially 
and non-axially symmetric vibrations of thick-walled 
hollow poroelastic cylinders and wave propagation in 
an empty cylindrical bore are investigated (Malla 
Reddy and Tajuddin, 1999, 2000) in the frame work 
of Biot's theory of wave propagation in fluid saturated 
porous media. Solorza and Sahay (2004) developed 
expressions for torsional resonance and temporal 
attenuation frequencies for fully saturated poroelastic 
circular cylinders. Poromechanic analysis of a fully 
saturated transversely isotropic hollow cylinder is 
made (Kanj et al, 2003; Kanj and Abousleiman, 
2004). Ahmed Shah (2008) studied axially symmetric 
vibrations of fluid-filled poroelastic circular cylindrical 
shells. In Ahmed Shah (2008), the frequency equations 
of axially symmetric vibrations propagating in a fluid- 
filled and an empty isotropic poroelastic bore, are each 
investigated for pervious and impervious surfaces. 
Radial vibrations of thick-walled hollow isotropic por- 
oelastic cylinders are studied (Tajuddin and Shah, 
2010). Flexural wave propagation in coated 
poroelastic cylinders with reference to fretting fatigue 
is studied (Ahmed Shah, 2011). In Tajuddin and Shah 
(2010), the dissipative coefficient is neglected due to the 
mathematical complexity. To the best of the authors’ 
knowledge, radial vibrations in thick-walled trans- 
versely isotropic cylindrical bone in the presence of dis- 
sipation are not yet investigated. In the present 
paper, bone with marrow is modeled as a transversely 
isotropic poroelastic cylinder filled with fluid. The per- 
tinent governing equations in the case of radial vibra- 
tions are derived. The frequency equations in the 
presence of dissipation are derived in the cases of 
bone without marrow and bone with marrow. In a 
limiting case, the attenuation as a function of the 
ratio of radii is computed for various anisotropic 
ratios. 

The rest of the paper is organized as follows. In 
Section 2, cylindrical bone without marrow is con- 
sidered. In Section 3, first behavior of waves in 
marrow independently from the solid bone is investi- 
gated, next the bone with marrow is examined. 
The limiting cases and numerical results are 
discussed in Section 4. Finally, conclusions are given 
in Section 5. 


2. Cylindrical bone without marrow 


Consider cylindrical coordinate system (r,0,z) with z- 
axis along the axis of poroelastic cylindrical bone. Let 
the bone be homogeneous and transversely isotropic. 
Assume that z-axis 1s along the axis of rotational sym- 
metry. Hence, in this case eight constants are involved. 

The constitutive stress-strain relations for a trans- 
versely isotropic poroelastic solid (Abousleiman and 
Cui, 1998) are given below. The notations are adopted 
from the said reference. 


Oy Mi Mp M3 0 0 0 Cpr 
O6 My Mi M3 0 0 0 eo 
oz | | Mi3 Mi M3 0 0 0 zz 
tw} | 0 0 0 M4 0 0 İlya 
Toz 0 0 0 0 Ms 0 Yoz 
t 0 0 0 0 9 Mall 

o 

o 

y 

= 0 P» 
0 
0 
p= Mle— oley + ego) — t ez] (1) 


In equation (1), ej and y; are strain components; oj 
and t, are stress components in the cylindrical coord- 
inate system, respectively (Abousleiman and Cui, 1998); 
pis the pore pressure, o and a’ are Biot's effective stress 
coefficients in the isotropic plane (r-0 plane) and in the 
z- direction, respectively; M is Biot’s modulus, € is the 
variation of fluid content per unit reference volume, 
and Mj, are components of the drained elastic modulus 
which depend on E, E, v v, G and G'. E and v are 
drained Young's modulus and Poisson's ratio in the 
isotropic plane, £ and v' are similar quantities as that 
of E and v pertaining to the direction of the axis of 
symmetry,G and G' are the shear modulus related to 
the direction of the isotropic plane and axis of sym- 
metry, respectively. The symbols œ and a’ can be 
expressed in terms of Mj, and the bulk modulus of 
solid constituents K, and computed as in 
Abousleiman and Cui (1998). Young’s modulus (£) 
and Poisson’s ratio (v) (which leads to all necessary 
Mj), M and K, are suffice to compute all the above 
parameters. For given anisotropic ratios of Np <£ 
and N, =X (Kanj and Abousleiman, 2004), E' and v 
can be determined. Different Ng and N, ratios define 
different degrees of anisotropy. 

The equations of motion in transversely isotropic 
poroelastic solids in the presence of dissipation (5) are 
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not readily available in the literature, hence the same 
are derived from the constitutive relations 
(Abousleiman and Cui, 1998) and equations of equilib- 
rium, and are given below: 

Ma 9 / Ms? 


Mii — Mio a 19 1 
+ + 5 5 
2 02 rər r? 2r? 09? 2 dz? 
Mi, —M 0 1 
+ Mp + Mo? + = E d Ur 
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Ges) rt 


2 drdz or 
= ri + U y + pte = 
— 9? pur + pU, at Ur 
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rd r 00 (5 T ;) "Uo p 08 r 08: 
9 a 
= ag 2e + 23 Us) — ba e L9), 
P oM a (5 " Ju aM up Mel Pu- 
or r r 000z dz? 


0 
= ag onn: + p22Uz) — bare — U.), 
(2) 


where à (u,, ug, Uz) and U (U,, Us, Uz) are solid and fluid 
displacements, e and e are the dilatations of solid and 
fluid, respectively; oj, are mass coefficients. 

Consider the transversely isotropic thick-walled 
cylindrical bone with inner radius r;, and outer radius 
r2. Non-zero displacement components of solid and 
fluid for radial vibrations are u(u,,0,0) and 
U (U,, 0,0), respectively. These displacements are func- 
tions of r and time ¢ only. Equations (2) in cylindrical 
polar coordinates when the solid and fluid displacement 
components u, and U, are independent of 6 and z 
reduce to 


1 de de 
mm (v - =) ur + (Mir + Mo? + Misi) aM 


a ð 
= ag onm Tp» U,) = be (ur = U,), 


n: 
or or 


2 
= lor + pU) -b u-u) O) 

In equation (3), V? is the Laplacian. The solid and 
fluid displacement components can be evaluated from 
equation (3) representing plane harmonic vibrations. In 
this case, the displacements of solid u, and fluid U, in 
the radial direction are taken as follows: 


ur (r,t) =f — (4) 
U,(r,t) = Frye 


Substituting equation (4) in equation (3), then equa- 
tion (3) becomes 


(Mii + Me*)Af —aMAF = —w(Ky f + KP), 
—aM Af + MAF = (Kaf + KnF) 


(5) 
Solving the above equations, there is obtained 
F=XjAf -Xf (6) 
where 
2 —MMıı 2 MKıı -- aM Ki» 
l X o*MMmj; — aM») © MKp -- 4MM»' 
& id 1 
A—— 4———— (7) 


d? rd r 
Substituting F into equation (5), there is obtained 


Af Of +Xif =0 (8) 
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where 


xig? (Mii + Mo2)M» + MM), - &MKı3 
$ MM > 


Ki Ky) — Kj 
3829 2). 
que ( MM, 


Equation (8) can also be written as 


ALA HEDA (7) 0 (9) 


where & and & are given by & + & = X2, && = Xj. 
Equation (9) can also written as 
(A+é)f=0, 


(A+8)f=0 (10) 


Solving equation (10), there is obtained f for & and éz. 
Similarly, F can be obtained. Substituting these in equa- 
tion (4), the displacement components are as follows. 


u(r, t) = (cii (ir) ke Yılın) + e3Ji (Ear) + c4 Yi lére, 
U,(r, t) = —(c1 8; (ir) + 381 Yi (ir) 


+ c385J 1 (Ear) + c405 Yi (£r))e'^' (11) 


where cı, c5», c and c4 are arbitrary constants, w is the 
frequency of wave, J, and Y, are Bessel functions of 
first and second kind of order n, respectively, and 
0j, £j (j = 1,2) are 


g 
5 = (ons #aMKı3)— 5 MMu 


x(MKi;--auMK»), (j=1,2), 
|a + L3? (4-7 L3)? 
im La ° 7 La Í 


Li = o*((Mii + Mo))K» + MKi;  20MK)), 
L; = (w*((Mii + Me?) Kaz + MK şi + 24M Kio) 


1 
= 4o MM (Ki Ko = Kj», 


Ly = MM1) (12) 
In equation (12), Ki, = p11 — PA Kio = p12 +4 and 
Ky = py — Ë. 

Using these displacement components into stress- 
displacement relations (Biot, 1956), the relevant stress 
oy, and fluid pressure p are obtained, which are 


Ow = (1 A11 (F) + eii) + es Ais(r) + caAua(r))e"", 
p = (An) + @An(r) + 6342(r) + c4Ara(r))e 
(13) 


where 
2 3 
duit) ed — a) 


-iMi t Mio 4-2Maà(a — 1) J (ér) 
-F(Mói(a — 1) - CM + Mala — 1))&i Ji (5r), 
Azı) ME) cM) iir) MR c M Eu 


Aj2(r), A»»(r) are similar expressions as Aj1(7), Azi(r) 
with J1, J2 replacedby Y1, Y», respectively, 
A13(r), 473(r) are similar expressions as A 11(r), Azi(r) 


with £, 6; replacedby &>, 6>, respectively, 


Aj4(r), A24(r) are similar expressions as A 11(r), Azi(r) 


with J;, J5,£1,6; replacedby Y,, Y5,£»,ó», respectively 
(14) 


The stress-free boundary conditions in the case of a 
permeable boundary of thick-walled hollow (without 
marrow) bone are 


Oy tp=0, p=0 atr=r,andr=r (15) 


In the case of impermeable boundary, there follows 


9 
Or -- p — 0, = =0 atr—r,andr—r> 


(16) 

Eguations (13) and (15) together give a system of 
four homogeneous equations in four constants 
C1, C2, C3 and c4. In order to obtain a non-trivial solution 
of this system, the coefficients matrix must be singular. 
Accordingly, the frequency equation for permeable 
boundary is 


Ayı) Air) Aii) Air) 
Ax(ri) Anr) Az(rı)  Azalrı) 20 (17) 
Air) Anr) Ais(r2)  Aia(r2) 
Azı) A»xXro) Az3(r2)  Ao4(r2) 


Similarly, the frequency equation for impermeable 
boundary is 


Au(ri) Arı) Arr) Au) 
Bari) Barı) Bz(rı) Bon) cs (18) 
Ayır) Ai) Ais(r) Ayla) 
By(r;) Boro) Bos(r2) Bul) 
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where, 


2 
By (r) = - (M8 + aM) & Jx(&ir), 


B»(r) are similar expressions as B»j(r) with J 
replaced by Y», respectively, 

B>3(r) are similar expressions as Bo\(r) with &, ôi 
replaced by £», à», respectively, 

B»4(r) are similar expressions as B>1(r) with Jo, &, 6) 
replaced by Y, &,52, respectively. 

If the transverse isotropy effect is neglected, that is, 
poroelastic cylindrical bone is isotropic, the frequency 
equations (17) and (18) coincide with that of the paper 
(Tajuddin and Shah, 2010). 


3. Cylindrical bone with marrow 


Bone marrow is a key component of the circulatory 
system, producing the blood cells that support the 
body’s immune system. It is a kind of fluid resembling 
blood, which behaves like a non-Newtonian thinning 
fluid. For the radial vibrations in the marrow, the dis- 
placement potential function @ satisfies the wave 
equation: 


rp lap 199 
or rər c 02 


(19) 


where cy is the velocity of the sound in the marrow. 


The solution $ is given by ¢ = Jo(ro/ cy )e'^'. The pres- 
sure p(— pr) and radial displacement D,(— ©) are 
computed. The ratio of the pressure to radial displace- 


ment is 


Df _ pro» Jo(rof cy) 
D, — ofegJi(rofcy) (20) 


The boundary conditions to be satisfied on the sur- 
faces r = rı and r = r to be stress-free in a bone filled 
with marrow for a permeable boundary are 


WEP P, poo atr=nandon +p =0. 
u D, 
p=0atr=r 9 


The first boundary condition presents a condition at 
the fluid-solid interface, while the other conditions 
remain the same as in the earlier section. The boundary 
conditions for the impermeable boundary are 


rr 0, 
2 um Pf v=, atr = rı ando,,+p=0, 
u D; or 
d 
P =0atr=r (22) 
or 


For a non-trivial solution, the matrix of coefficients 
must be singular. Accordingly, the frequency equations 
in the case of a permeable boundary is 


Airy) Ap) Aig) Aja) 
Adri) Ayrı) Azri) AG) E (23) 
Air) Anr) Alr)  Aj(r2) 
Ax (r2) Anr) Aos(r2) Ax (r2) 
where 


An) = -Anw lro/ ce) + p Jo(rof cy) Ei), 
Aylı) = —Ara(r)o/ci(ro/ cp) + p Jo(roof cy) Yi (ir), 
Aylı) = -Arw cloc) + p Jo(roof cr) (Eor), 
Ay) = —Au(r)o/ejfi(rofcy) + p Jo(rof cy) Yi (E20), 
Abr) = Ar), j—1,2,3,4 

(24) 


The frequency equation in the case of an imperme- 
able boundary is 


Ai) A'sri) A'gri) A'a(ri) 
Ba(ri) Bo(ri Barı) Bon) = (25) 
Air) Anr) Ail) A42) 
By (12) Boro) Bor) Ba(r2) 


The frequency equations (17), (18), (23) and equation 
(25) will be examined numerically in a limiting case. 


4. Numerical results 


In the presence of dissipation (b), each element of fre- 
quency equations (17), (18), (23) and equation (25) con- 
sist of Bessel's functions and the argument of Bessel 
functions are complex valued, hence giving mathemat- 
ical complexity in numerical evaluation. To reduce 
complexity, the limiting case is considered. The limiting 
case when Ł « 1, (i.e. &ri, &ro, £j, £yr; >> 1) is con- 
sidered so that asymptotic approximations for Bessel's 
function (Abramowitz and Stegun, 1965) can be used. 
Under this condition, cylindrical bone reduces to thin 
cylindrical poroelastic shell type bone. The frequency 
equation in the case of bone without marrow for a per- 
meable boundary is 


Cu(rn) Cpo(n) Cin) Cut) 
Cari) Co(n) Cm») Cani) 20 (26) 
Ci) Ci(r2 Cis) Cur) 
Cr(r2) Cl) Cr) Ca(r2) 
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where Ci) = cpr) sidir), f= 1,2; 121,2, 
3,4, i— J—1. The expressions for cj(r), cir), 
J=1,2,7=1,2,3,4, are expressed in Appendix A. 
Similarly, the frequency equation in the case of bone 
without marrow for an impermeable boundary is 


Cu(u) Crt) Ciao) Cun) 
Da(r) Do»(ri Dori) Do4(ri) e en 
Cu(r) Ciro) Cız) Cir) 
Do(ro) Do»(ro) Dos(ro) Da(r2) 


where Dor) = dar) + id, (r), 121,2, 3, 4, i= J—1. 
The expressions for dor), d)(r), L= 1, 2, 3, 4, are 
expressed in Appendix A. Under the condition, 
Ae 1, cylindrical bone with marrow reduces to thin 
cylindrical poroelastic shell type bone with marrow. 
Skull, nasal, hip bone, shoulder blades and ribs are of 
shell type. These bones are thin and generally curved 
and have red bone marrow rather than both red and 
yellow bone marrow. The frequency equation (23) in 
the case of bone with marrow for permeable bound- 
ary after substituting the asymptotic approximations 
would be 


Piri) Pit) Pzr) Pun) 
Py(r) Pari) Por) Poi) 20 (28) 
Cu(r) Ciro) Cız) Cir) 
Ca(r2) Coxro) Caz) Cali) 


where Pır) = pir) + ipy(r), Por) = Cor), 1— 1, 2, 
3,4,i—N—1. The expressions for pur). pr). 
l= 1,2,3,4 are given in Appendix B. Similarly, the 
frequency equation for an impermeable boundary is 


Pyu(r Piri) Pisin) Pur) 
Da(ri)) D(ri) Dzi) Doa(nri) - (29) 
Cu(r  Cio(2) Cir) Ciro) 
Do(ro) Do(ro) Dos(ro) Doa(r2) 


where the expressions Pı/(), D2/(r), Cyi(r), | = 1,2, 3,4 
are given in Appendix A and B. The frequency equa- 
tions (26)-(29) are investigated numerically. The mater- 
ial constants M11, Mi», M13 involves E, E', v and v’. The 
values of Young's modulus (£), Poisson's ratio (v), 
Skempton pore pressure coefficient (B), undrained 
Poisson's ratio (v,) and dissipative coefficient are 
respectively taken to be 14.58 Gpa, 0.32, 0.4, 0.33, 
and 0.17 x 1076 Gpa/m/, respectively, as suggested in 
Cowin (1999). In Steck et al. (2003), it is assumed that 
the anisotropic ratios of bone vary in the neighbour- 
hood of one. The bulk modulus of solid constituents 
(K,) and the Biots modulus (M) can be determined 


from the following relations (Kanj and Abousleiman, 
2004): 


pe BA + v) E 
* 8B(1 — 2vY(1 + vu) — 9(v, — v)’ 
B(1 — wW)(1 + v) E 

9(1 — 2v,)(1 + vv, — v) 


The values of mass coefficients pj, are computed as 
in Nowinski and Davis (1972). 

The velocity of sound in marrow (cy) and density of 
marrow (oy) are taken to be 93m/s and 900 kg/m?, 
respectively (Malla Reddy et al., 2011). Employing 
these values in the frequency equations, an implicit rela- 
tion between complex frequency (w) and the ratio 
(g — 2) is obtained. The attenuation coefficient (Q^!) 
(Solorza and Sahay, 2004) is given by 


M= 


_ 2Im(o) 


=i 
9 = Re(o) 


(30) 


The phase velocity (cj) and group velocity (c,) are 
given by 


cy = Real part (w/k) (31) 
dw 
C= dk (32) 


Attenuation against the ratio (g) is computed in the 
cases of bone without marrow and bone filled with 
marrow for permeable and impermeable boundaries. 
Phase velocity and group velocity are computed as a 
function of wavenumber in the case of bone without 
marrow for a permeable boundary. The numerical 
values are depicted in Figures 1-6. The notations 
Bwm, Bwom, Per and Imp used in the figures to repre- 
sent the bone with marrow case, bone without marrow 
case, permeable boundary case, and impermeable 
boundary case, respectively. Figure 1 depicts the vari- 
ation of attenuation against the ratio (g) when the 
anisotropic ratios Ng = land N, = 0.5 in the cases of 
bone without marrow and bone with marrow for per- 
meable and impermeable boundaries. From this figure, 
it is clear that the attenuation values in the case of bone 
with marrow are, in general, greater than that of bone 
without marrow in the case of a permeable boundary. 
Also, it is seen that the attenuation values in the case of 
bone with marrow are, in general, less than that of bone 
without marrow in the case of an impermeable bound- 
ary. That 1s, the trend is reversed for an impermeable 
boundary. Figure 2 shows the attenuation values 
against the ratio (g) when the anisotropic ratios 
N,=1 and Ng —0.5 in the cases of bone without 
marrow and bone with marrow for permeable and 
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Attenuation 


—— Bwom, Per 
—A— Bwm, Per 
—4— Bwom, Imp 
—*— Bwm, Imp 


Ratio (g) 


Figure |. Variation of attenuation against the ratio (g) when the anisotropic ratios Ne = | and N, = 0.5 in the cases of bone 
without marrow and bone with marrow for permeable and impermeable boundaries. 


12 — —e— Bwom, Per 
—*— Bwm, Per 
10 + —A— Bwom, Imp 
—— Bwm, Imp 
> 8 
E 
3 
E 6 
* 4 
2 
0 
1.1 1.5 
Ratio (g) 
Figure 2. Variation of attenuation against the ratio (g) when the anisotropic ratios N, = | and Ne = 0.5 in the cases of bone 


without marrow and bone with marrow for permeable and impermeable boundaries. 


impermeable boundaries. From Figure 2, it is observed 
that attenuation values in the case of bone with marrow 
are, in general, less than that of bone without marrow 
in the case of a permeable boundary. It is also seen that 
attenuation values in the case of bone with marrow are, 
in general, greater than that of bone without marrow in 
the case of an impermeable boundary. From Figures 1 
and 2, it is clear that the nature of surface does have an 
influence over attenuation values in the cases of bone 
without marrow and bone with marrow. Figures 3 and 
4 depict the variation of phase velocity and group vel- 
ocity with wavenumber (K) in the case of bone without 
marrow for a permeable boundary when the ratio 
g — l.l. From Figures 3 and 4, it is observed that 
phase velocity and group velocity values at anisotropic 
ratio N, — 0.5 are, in general, greater than that of 
values at N,=0.75 when Ng—l and values at 
Ng — 0.5 are, in general, less than that of values at 


Ng — 0.75 when M, = 1. Figures 5 and 6 depict the 
variation of phase velocity and group velocity with 
wavenumber (k) in the case of bone with marrow for 
a permeable boundary when the ratio g — 1.1. From 
Figures 5 and 6, it is seen that phase velocity values 
at anisotropic ratio N, — 0.5 are, in general, greater 
than that of values at N, = 0.75 when Ng—1 but 
group velocity values at anisotropic ratio N, — 0.5 
are, in general, less than that of values at N, — 0.75 
when Nç = 1. Similar conclusions are made in the 
case of anisotropic ratio Ng = 0.5, 0.75 when N, = 1. 
From Figures 3-6, it is clear that, for given anisotropic 
ratios, the trend of phase velocity and group velocity is 
the same in the bone without marrow case but different 
in the bone with marrow case. From all the figures, it 1s 
found that dispersive behaviors in bone without 
marrow and bone with marrow cases are different in 
phenomena. 
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Figure 3. Variation of phase velocity with wavenumber (k) in the case of bone without marrow for permeable boundary when the 
ratio g = l.l. 
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Figure 4. Variation of group velocity with wavenumber (k) in the case of bone without marrow for permeable boundary when the 
ratio g = l.l. 
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Figure 5. Variation of phase velocity with wavenumber (k) in the case of bone with marrow for permeable boundary when the ratio 
g-l.. 
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Figure 6. Variation of group velocity with wavenumber (k) in the case of bone with marrow for permeable boundary when the ratio 


g- l.l. 


This discrepancy can be exploited in the medical 
diagnosis to examine the level of bone. When the 
ratio of outer and inner radii of bone is fixed, the 
phase velocity values are, in general, greater than that 
of group velocity in both the cases of a permeable 
boundary at various anisotropic ratios. Attenuation, 
phase velocity, and group velocity values are affected 
with anisotropic ratios in either case. 


5. Conclusion 


Thick-walled cylindrical bone filled with marrow is 
modeled as a transversely isotropic thick-walled cylin- 
drical poroelastic cylinder filled with fluid in the frame- 
work of Biot’s poroelasticity theory. The pertinent 
equations of motion are derived. A comparative study 
has been made between the bone without marrow and 
the bone with marrow for the cases of permeable and 
impermeable boundaries. The limiting case, when the 
ratio of thickness to inner radius is very small, is dis- 
cussed. In this limiting case, thick-walled cylindrical 
bone reduces to thin shell type bone. The available 
data of bony materials has been used for numerical 
evaluation. Attenuation versus ratio of outer and 
inner radii of bone at various anisotropic ratios are 
computed. Phase velocity and group velocity are com- 
puted as a function of wavenumber at various aniso- 
tropic ratios in the cases of bone without marrow and 
bone with marrow for a permeable boundary at various 
anisotropic ratios. From the figures, it is observed that 
the nature of the surface does have an influence over 
attenuation values in the cases of bone without marrow 
and bone with marrow. From the figures, it is also clear 
that, for given anisotropic ratios, the trend of phase 
velocity and group velocity is the same in the bone 
without marrow case but different in the bone with 


marrow case. From all the figures, it is found that dis- 
persive behaviors in bone without marrow and bone 
with marrow cases are different in phenomena. This 
discrepancy can be exploited in the medical diagnosis 
to examine the level of bone. When the ratio of outer 
and inner radii of bone is fixed, the phase velocity 
values are, in general, greater than that of group vel- 
ocity in both cases of permeable boundary at various 
anisotropic ratios. Attenuation, phase velocity, and 
group velocity values are affected with anisotropic 
ratios in either case. 
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Appendix A 


eu(r)- — (M(1 — aa HUS + 2Mo(a — 1) 


+ Mio) Dy)bi — (My + Ma(a — 1))a» 
+ M(1 — a)a3)bo, 


2 1 
eix(r) = 7 (M(1 — a)a4 ; (Mu -- 2Ma(a — 1) 


+ Mio) D;)bi — (My, + Ma(a — 1))a» 
— MÜ — o)ac)b», 


—2 

ca (r) = , Ma — aMD,)b, + (aMap + Maa)b», 
—2 

cox(r) = — Mas + aMD>)b, + (aMas + Mag)bo, 


c13(r), cı4(r) are similar expressions as c1;(r), cı>(r) 
with ai, o, a3, bi, b2, Dı, Do replaced by a7, dg, a9, 
b3, b4, D3, Da, respectively, 


a = ài (Tii cos(7/8) + T1» sin(z/8)) + 013(T15 cos(z/8) 
— Ty sin(z/8)), 
a = Li(cos Ti(T13 cos(z/8) + Ta sin(z/8)) 
— sin Ti(T14 cos(7:/8) — T3 sin(7:/8))), 
az = (T\3 cos(z/8) + Tıg sin(1r/8)X(611/, cos T; 
+ àol, sin T1) — (T13 cos(7:/8) + Ta sin(z7/8)) 
x (811lj cos Tj + 02/1 sin T1), 


a4 = ó11(T21 cos(zt/8) + Tz sin(z1/8)) + 512(T24 cos(7/8) 
— T» sin(z/$)), 
as = lh (cos Ti(T53 cos(7/8) + T54 sin(77/8)) 
— sin T3(T54 cos(7:/8) — T>3 sin(7:/8))), 
à = (Tas cos(z/8) + Taa sin(z/8)) 
x (8114, cos T, + 0joli sin Ti) — (T54 cos(z/8) 
— To sin(z/8))(61,/) sin Tı — 613/1 cos T3), 


b, = | -cos'Gr/8) — sin? (1/8)! sin(377/4), 
linr 


bo = ET — sin^(x/8)) !, 
linr 


Dı = Tii(cos(7/8) + Ti sin(zr/8)), 
D» = T»(cos(7/8) + To sin(z/8)), 
h — (di di, T, 5 z Tan (2). 


|. (Mii + Mo?) px + Mpi + 20M pio 


di b cos T: 
XMa MM yp 19 me» 
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FTN —ob((Mi + Mo?) + M(14- o) 
27 (Mı + M9)M — a M3) 


1 d 
b — (d 4- df, T; = Tan" (5 ], 
2 d 


+h sin T5, 


ds = (Mii + Me?)p> + Mpii + 2aMpı3) 
b2 
tO Mo? + M —2aMy) 
— 4o (piv 022 — i5), 
2b 2 
d4 = Aob(pii + 2012 + p2)M 1 — p Wn + Ma^) 


+ Mpii + 2aMpı3 Mı + Mo? + M — 20M), 
© Mpu +aMpı2 — MMF cos(2Ty))/o? 
ET M pio + & M p»» 
© MM uil sin(2T))/@ 
© Mp kaMp» 
Tı, = sin(/; cos Tir) cosh(/; sin Tir 


E) 


— 


— cost, cos Tır) cosh(/, sin Tir), 


E 


= cos(/; cos Tir) sinh(/; sin Tr) 

+ sin(/, cos Tir) sinh(/; sin Tir), 
T3 = sin(/; cos Tir) cosh(/; sin Tir) 
+ cos(/ cos Tır) cosh(/; sin Tir), 


3 
E 


= cos(/; cos Tr) sinh(/, sin Tir 


— 


— sin(/, cos Tr) sinh(/, sin Tir), 


15 
Thy = x1 + X3 + SI (Cy1x1 — y2x2) — Gayi + y2x4)), 


Tx = X4 X2 


15 
Sr ((y1xa — y223) + Orayı + y2x1)), 


L5 
T33 = x3 — x1 + g OS + yoxa) — (X11 + ¥2X2)), 


15 
To4 = x4 + X2 + 8r (xa — 2x3) — Gayo + yoxi)). 


xı = cos(/, cos Tır) cosh(/; sin Tir), 
xa = sin(/; cos Tir) sinh(/; sin Tir), 
x3 = sin(/, cos Tir) cosh(/, sin Tr), 


x4 = cos(/; cos Tır) sinh(/; sin Tır), 


yı = cos Ti, y2 = sin Tj, 


T31, T32, T33, T34, Ta, Ty, T43, Taa are similar expres- 
sions as Tu, Ti», Tia, Jas T», Tr, T3, T4 with lh, Ti 
replaced by 4,73, respectively, a7,ag,a9 are similar 
expressions as d1,405,4d3 with Tu, Ti», Tis, T\4, li; 
Tı,bı,b2 replaced by T31, T32, T33, T34, l3, T3, 53, ba, 
respectively, 


b3,b4 are similar expressions as bj,b5 with h 
replaced by h, 


D3 = T3 cos 7/8 + T32 sin 7/8, 

Dı = Ty cos 1/8 + Ty sin 7/8, 
419, 411, 412 are similar expressions as a4, ds, aç with 
Tai, Tx, T23, To4, h, Ti, bi, b2 replaced by Ty, T», 
T43, T44, l5, T3, b3, b4, respectively, 521, 022 are similar 
expressions as 611, ô12 with /,, 7, replaced by h, T3, 
respectively, 


a3 = ój (T5 cos 7/8 — T, sin 7/8) 
— ôi2(T11 cos 2/8 + Ti» sin 77/8), 
daya = L»(sin Tı(Tı3 cos x/8 + Tıg sin 7/8) 
+ cos Ti(T314 cos 2/8 — T3 sin 7/8)), 
ais = (Tia cos z/8 — Tız sin z/8)(814 L4 cos Tj 
+ 015 L4 sin T1) + (Ti cos 2/8 + Ti4 sin 7/8) 
x (6, Lı sin Tj — 124 cos T3), 
Ds = Ti» cosz/8 — T sin 7/8, 
D; = T35 cosz/8 — T3 sin 7/8, 


Cızlr), C4(r) are similar expressions as cjj(r), cio(r) with 
41,02,03,D1,b5, Dı, D2,ô11,ô12 replaced by a7, ag, do, 
b3, b4, D3, D4, 521, 522, respectively, cos(r), co4(r) are 
similar expressions as coj(r),co(r) with  a4,as, a6, 
bi, b2, Di, D2, 511,512 replaced by aio,a11, a12,53, ba, 
D3, D4, 521,522, respectively, 


1 d, 
(Ed, Ts tan (2). 


| (Mii + Ma’) p22 + Mpii + 20M pi2) 


d. L T: 

5 2Mii 2COS 12, 
— M Mo?)4- M +a? M? 

d= cb((M; + Mo?) + M +a M?) bsin To, 


dy) AMA oğan 

ahı +2Ma(a—1)+ Mi2) D3)bi 

— (M1 + Ma(a İla Mİ aan bs 
cir) -4M0 — ajay +n +2Ma(a— 1) 


+ Mi)De)bi (My, + Mala — 1))ai; 
— M(1 — a)ais)bz, 


—2 
eur) -- Mas -FaMDs)b; -- (aMaya + Maıs)bo, 


2 
cr) — (Marg +aMD,)bh; + (a Maj; + Majg)ho, 
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419, 429, 42, are similar expressions as 13,444,445 with 
Tiu, Tio, Tis, Ti4, hi, Ti, 811,812, Ds replaced by T31, T32, 
T33, T34, 5, T3, 521, 522, D7, respectively, a22, 423, 424 are 
similar expressions as 446,4175,d]g with Tzn, Tx, 
T53, T24, lh, T1,811, 612, De replaced by 141, T», Tas, 
TAa, l3, T3, 621,622, Dg, respectively, 


ayı = l(cos Ti(T53 cos 7/8 + T54 sin 7/8) 
— sin Ti(T54 cos x/8 — Tzs sin x/8)), 
aig = (T54 cos 1/8 — T>3 sin 7/8) 
x (5114, cos Tı + 124 sin Ti) 
+ (T53 cos 2/8 + Tag sin 7/8) 
x (5114; sin Tı — 6124) cos T3), 
Do = Ti4 cos z/8 — Tız sin 7/8, 
De = T34cos 27/8 — T33 sin 7/8, 


c3(r), c14(r) are similar expressions as c,(r), ci>(r) with 
413,014, 415, b1, b2, Ds, Do, 611,512 replaced by ayo, 
250, 4215 b3, b4, D}, Ds, 651 ‘i 055, respectively, Cyr), cor) 
are similar expressions as chi(r),ch (r) with a6, 


417,018, b1, b2, Ds, Dg, ô11,ô12 replaced by — a25, 493, 
dra, b3, ba, D7, Dg, 521, 62 respectively, 
2 : 
dy (r) = z ae + a M)(Dob cos T» — Dioh sin T5) 
+ M6p(Dıgh COS T» + Dol sin T>)), 
2 , 
dx(r) = z Mon -- aM)(Diil» cos T» — Diob sin T5) 


+ Móio(Diilo sin T» + Diolo cos T>)), 
2 : 
dr) = z (Môi + aM)(Dioh COS Tə + Dols sin Tə) 


= M6p(Dyh cos T^ = Dioh sin T2)), 


2. ; 
dr) = 7 (Mà +aM)(Di,h sin T» + Diol cos T2) 


— Móyi(Diib cos T» — Diob sin T>)), 


Do = Tx, cosz/8 + T»; sin 7/8, 
o = Tx) cos 1/8 — T», sin 7/8, 
11 = Ta cos 7/8 + Tag sin 7/8, 
2 = Tx cos 7/8 + Toz sin 7/8, 
3 = Ty cos 7/8 + Ty sin 7/8, 
4 = Ty cos z/8 — T4; sin 77/8, 
5 = Ty cos 7/8 + T44 sin 7/8, 
6 = T44 COS 7/8 + Ta sin 7/8, 


ee Miele 


dy3(r), dalr), b(r), d,,(r) are similar expressions as 
dın), der), di (r) d», (r) with h, T», Do, Dio, Dii, 
Di2, ô11,ô12 replaced by o 5, T5, Dis, Dia, Dis, 
Dıs, 521,522, respectively, 


Appendix B 


= ro ræ 
pu(r) 2 — Jı 2) cur) + per Ea Dı, 
Cf Cf Cf 
—eo fr ro 
pir) = —4 (=) cnlr) + se n Ea D», 
Cf Cf Cf 
—o , (ray , ro 
py”) = a (=) cr) + ojo "n (ns. 
Cf a 
/ p rw / 2 
pir) = =h — | e(r) + pw H- = ? Jb Do, 
Cf T. Cf 


pia(r), pia (r) are similar expressions as pıı(r), pio(r) with 
cu(r),ci(r), Di, D2,b1, replaced by — eis(n) c14(r), 
D3, D4, bs, respectively, 

Pia (D. pi4(r) are similar expressions as pjj(r), pi;(r) 
with cj(r), cia(r), Ds, Ds, b replaced by cr), c4), 
D7, Dg, bs respectively. 
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